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SECTION  1 


INTRODUCTION  AND  OVERVIEW 

The  objective  of  this  research  is  to  evaluate  the  effects  of  viscous 
and  tirbmence  phenomena  on  the  flowfield  behind  a  blast  wave  which  passes 
over  a  structure.  Viscous  effects  have  not  been  previously  considered  in 
evaluating  such  phenomena.  However,  ic  is  well  known  that  viscosity  has 
global  effects  on  steady  flows,  such  as  being  the  mechanism  that  leads  to  flow 
separation.  Therefore,  it  is  quite  likely  that  the  predicted  blast  wave 
loading  on  a  structure  is  affected  by  viscous  forces  in  the  flowfield  and  that 
if  is  important  to  include  these  effects  in  comparisons  with  laboratory  data. 
Moreover,  there  is  a  need  for  flowfield  predictions  that  resolve  boundary 
and  shear  layers  because  other  effects,  such  as  dust  lofting  and  transport, 
a re  influenced  by  these  details  of  the  flow. 

It  is  useful  to  divide  the  interaction  time  into  two  primary  regimes, 
aarly  and  late  time.  Wo  will  use  early  time  to  refer  to  the  interval  during 
whi'-n  and  immediately  after  the  blast  wave  passes  over  the  structure.  The 
late  time  period  corresponds  to  the  tim?  when  the  shock  front  is  far  from  the 
ibjoot  (compared  to  typical  dimensions  of  the  structure). 

The  highest  impulsive  loads  occur  during  the  early  time  regime.  The 
vi  irons  layers  are  generally  quite  thin,  and  the  flowfield  is  primarily 
lovis.'id.  However,  depending  on  the  configuration  and  flow  conditions,  the 
vi  icons  layer  might  separate,  resulting  in  a  radically  different  flowfield 


structure,  us  illustrated  in  figure  1-1,  wh’  .h  shows  i  schematic  of  a  blast 
wave  flowfield  passing  over  a  vehicle.  The  asymmetric  pressure1  loads  have 
resultant  forces  and  moments  which  may  lamage  the  vehi  le,  move  it  laterally, 
or  turn  i‘  iver.  The  resulting  mot  i  in  nr  lamge  varies  depending  on  tne 
history  if  the  developing  t  low!  ield,  in.  i  idin;  whether  or  not  tne  flow 
separates . 

'lata  show  that  the  rate  at  which  a  Mowtield  run  change  from  attached 
f  i  separated  flow  scales  with  the  Sf.rotiha  1  number,  st  -  u  where  u  is  the 
■■lean  stream  ‘low  velocity,  -  is  a  typical  length  of  the  structure,  and  x  is 
fne  tone  interval.'-  For  impulsive  flowfield  changes,  St  -  0.17.  For  a  5  psi 
iverpressure  blast  wave,  with  10  ft,  this  gives  i  =  0( 10-3)  s.  For  a 
typi-.al  vehicle  under  these  conditions,  the  angular  rotition  rate  due  to  the 
large,  asymmetric  side  loads  is  roughly  10°/s.  Thus,  the  vehicle  response 

is  on  the  order  of  0.1  to  10  s,  and  it  is  clear  that  separation  can  occur 
;  i ’  ■  *  1 y  enough  to  affect  the  dynamics  of  forces  and  moments  on  vehicles  and 
.  ‘  r.ji.t.  UTeS  . 

during  the  late  time  period,  tne  flowfield  structure  changes  more 
;r.iduully.  Thick  viscous  layers  develop,  becoming  comparable  to  or  larger 
'nut'  t no  neight  of  typical  structures.  The  turbulence  flow  may  loft 


'vhoni.it.  i  c  of  blast  wave  flowfield  interacting  with  a 

field  veh  i  1  o 


1 


onsiderable  debris  (pebbles,  rocxs,  dirt,  etr.)  into  the  flowfield,  thereby 
■>t,erf  err  i  ng  with  communi  cat  1  nns  and  possibly  damaging  structural  components. 

Rocont  1  y ,  the  II. S.  Army  ballistic  Research  Laboratory  supported 
,>S'M'\.h  b,'  Sutler  and  Lerrujinst^  on  the  inviscid  flowfield  associated  with 
las'  wave -sf  met  a  re  interactions.  The  present,  study  extends  this  earlier 
r-  to  include  vi  scous./turbulence  effects.  We  have  used  the  same  kinds  of 
mite  difference  methods,  solution  procedures,  and  computer  codes  as  Kutler 
nd  rernqui st .  We  studied  first  the  simple  viscous  flow  in  a  shock  tube, 
m  ,;h  serves  to  validate  the  treatment  of  the  viscous  terms  in  the  equation  of 
of  ion.  these  results  compare  well  with  shock  tube  data  and  "exact"  boundary 
i y ' * r  cheer /  tor  this  case.  We  have  also  made  both  viscous  and  inviscid 
• '  •  ill*  i  .ms  tor  a  blast  wave  flow  over  a  free-standing  cylinder  and  for  a 
I.:  cylinder  resting  on  a  Mat  surface.  The  present  results  for  a 
fve-st  ind i n.)  cylinder  compare  quite  well  with  previous  results  and  data.  The 
iw  i  half  cylinder  plate  is  a  prototype  of  the  geometry  of  a  common 

i  -or  ion  M.q.,  building,  vehicle),  which  is  an  object  on  the  earth's 
f  •  <  it  i nr  hides  the  additional  feature  of  the  boundary  layer  development 
>  ,  • ar.*d  witn  the  incoming  flow,  which  appears  to  he  the  cause  of  an 

,•  r...i;-,  ,■>;).( r ,, t ,> r t  f  1  u w  region.  This  result  appears  to  he  a  unique  phenomena 
••  !  ’  ,  ’  ‘  >  be  :)f  practical  interest. 

’’  ,  report  is  tr-j  mi  /••.i  ;  n  four  parts .  Section  ?  describes  the 

;  j  ,  t*  i  )••  . ,  mi  a  1  i  jiritnn,  and  t;ne  boundary  conditions  used. 

•  i  •  •••  ,  ■■  ••  *  *•>..  *.•-,*  r,'- .  ni  f  he  i  r  results.  The  flow  behind  a  shocx 

i  ‘  •'  |-  •-  r’  ‘  •  ‘  i  i'-  it*  i  shock  and  a  cylinder  are  described  in 

•  ••  ••  •  ■  '  •  ■  d  .■.,•!>  i  ns  r  hi-  resu  I  f  s  of  our  ana  I  y  s  1  s  of  the 

•  -  i  *.••••  i  .►>  .•  ►  w  <  v  I'd  i  ”  n  ;•>  t  m  toe  ground,  finally. 


SECTION  2 


ANALYSIS 


The  physical  problem  of  interest  is  the  interaction  of  blast  waves  with 
objects  protruding  from  a  surface  as  shown  in  Figure  1-1.  However,  as  a  shock 
moves  along  a  surface,  a  boundary  layer  develops  behind  it  and  shear  forces 
change  the  flowfield.  Therefore,  to  adequately  describe  this  flow  the 
compressible  mean  Navier-Stok.es  equations  in  two  dimensions  (2-D)  are 
employed.  In  this  study,  relatively  weak  blast  waves  are  considered  and  the 
him*"  frame  of  the  calculation  is  relatively  short.  Therefore,  the  blast  wave 
■s  approximated  as  a  constant  speed  shock  wave.  In  addition,  the  yas  is 
!,,)’"■!  to  be  thermally  and  calorically  ideal. 

GOVERNING  EQUATIONS 

rhe  Nav ier-Stokes  equations  in  Cartesian  coordinates  for 
t wo-d'inensional  flow  can  be  written  in  dimensionless,  conservative-law  form 
fir  a  perfect  gas  without  external  forces  as  follows:? 
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'  .  ■  ;  ’  • '  i  '  ,  •»  ,  -  i»  i  .t: .  . . .  i  s  shown  except  near  the  shock  . 

'•■•  •  :•  ;  .  t '  • f  v  i*  '".i  >.,.)•  ,t  i  in  ne  !  r  the  shore,  is  expected  because  the 

.  •  *!-:*•  ■ :  ‘  t  ,■  .  1 1  jo  r  1 1  mu  of  AdC/i)  (with  fourth  order  dissipation) 

•  i  . s  n,ir.ifteri  st.  i  cs  for  other  users .  Some  effort  has  been 

i  ‘  ‘  ’  i  m  n  r.  i  hi>  fnis  problem.  ^  One  preliminary  solution,  the  addition 

•  ,.•!  on  •  >r  |,.r  l  i  ;s  i  pat  i  on  near  the  shock,  has  been  successful  in  steady 

>"  i f  ■  ■  i!  liit  ions.  Mnf ortunately ,  this  modification  was  not  available  for 
*  .e'O'.en'  study. 

fnese  results  demonstrate  that  the  code  does  a  good  job  of  predicting 
i  dinar  boundary  layer  behind  a  moving  shock. 
i  ■  i r b i j  1  ent  Results 

[he  development  of  a  turoulent  boundary  layer  behind  a  moving  shock  has 
■tie)  ne.>n  studied  both  experimental  ly  and  ana  lyt  i  cal  ly .  Hartunian,  et  al.,lf) 
Measured  heat  transfer  in  shock  tubes  for  a  wide  range  of  shock  strengths. 
v:irtin‘;  made  measurements  at  only  two  pressure  ratios  (2.75  and  3)  and 
.  umpired  his  results  with  the  approximate  turbulent  boundary  layer  theory  of 
t  r- •  1  s .  1 

[he  turbulent  boundary  layer  studied  here  is  created  by  a  shock  with  a 
►•res  ;  ratio  of  2.75  moving  along  a  flat  wall,  as  was  the  case  for  the 
1  imi  e,ir  study .  The  pressure  ratio  of  2./5  is  somewhat  unique  in  that  the 
-  1  I’d  t  velocity  behind  the  shock  in  the  shock  frame  of  reference  is 

1  ;  m  i  •  i  ni,.  fluid  velocity  in  the  wall  (or  fixed)  frame  of  reference.  A 
n,‘  r>*  will  t -‘tripe ra f. u rt-  is  assumed,  which  is  confirmed  by  the  measurements 
Martin  and  Hartunian  for  shocks  of  this  strength .  Compared  to  the 
i  ai "a  r  *  low  c  i  leu ! at  ion,  much  smaller  time  steps  were  used  for  turbulent  flow 
lie*  there  i;  finer  mesh  spacing  af.  the  wall.  The  grid  spacing  at  the  wall 
it  mi-'c,  .  ,  r.>  j,, i  red  to  ohfiin  resu  1  ts  that  are  independent  of  further 


Laminar  velocity  profiles  in  shock  frame  of  reference 
(Symbols  represent  calculated  values,  line  represents 
Mirels'  self  similar  solution) 
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a  function  of  nond i mens i ona 1  distance 
laminar  boundary  layer  behind  shock 


i . 1 . ?  Lam j nar  Resul t s 

Our  predictions  of  a  laminar  boundary  layer  are  compared  to  the 
"analytic"  results  of  Mirels  in  this  section.  We  will  follow  Mirels' 
transformation  of  the  unsteady  problem  of  a  shock  wave  moving  over  a  wall  to 
that  of  a  steady  flow  by  presenting  our  results  in  the  shock  fixed  frame  of 
reference.  Mirels  used  the  Prandtl  boundary  layer  equations  and  applied  a 
similarity  variable  which  is  a  function  of  the  freestream  velocity,  fluid 
parameters,  and  distance  behind  the  shock.  The  resulting  equation  is  an 
ordinary  differential  equation  which  Mirels  integrated  numerically. 

Martini  an,  et  al.,^(]  found  good  agreement  between  their  measured  heat  transfer 
rates  and  those  calculated  with  the  equations  derived  by  Mirels.  They  also 
showed  that  the  effects  of  density  and  viscosity  variations  are  small  for  Mach 
nuni>-,rs  less  than  2.  These  results  shew  Mirels'  theory  to  be  a  good  model  for 
laminar  flow. 

the  laminar  flow  calculations  hc  ve  a  shock  Mach  number  of  1.58.  The 
details  of  the  flow  calculation  are  given  in  Table  3-2.  After  60(1  steps  the 
shock  is  6.2m  from  the  leading  edge,  ard  all  the  results  we  will  show  are  for 
this  time.  Figure  3-3  shows  lines  of  constant  Mach  number,  and  qualitatively 
■narks  the  edge  of  the  boundary  layer.  It  is  easy  to  distinguish  the  two 
’'■'jinns  of  the  flow.  Figure  3-4  shows  good  agreement  between  the  calculated 
velocity  profile  with  that  predicted  b>  Mirels1  similarity  solutions.  In  this 
figure,  the  ordinate  is  Mirels1  similarity  parameter  and  the  velocities  are 
""'■i-.  i  red  r  -  *  5  a  t  i  ve  to  the  shock.  Thus  i  /ue  =  2  is  the  wall  and  u/ue  =  1  is  the 
‘  ?  '■"■im.  Additionally,  throughout  t  he  flow  the  edge  velocity  normal  to  the 

wi'l,  a  ■  t  n ,  )•  j  gh  not  shown,  has  the  correct,  behavior  and  is  directed  toward  the 
will  o-ar  the  shock  and  away  from  the  v.all  near  the  leading  edge.  Figure  3-5 
i  i  ■  o'lipari  son  of  the  computed  frictim  coefficient  with  that  predicted  by 


r  he  most,  suitable  tesf  tor  .1  -  .1 1  ulation  of  1  am  i  n  a  r  flow  behind  a  shock  is  a 


comparison  with  Mirels1  similarity  solution.1*  Mirels'  results  have  been 
compared  favorably  with  experimental  data,l()  and  the  analysis  required  to 
compare  to  the  work  1 qu  ,  t-j  si  ru  i.jht  forward  .  Comparison  with  experimental 
lata  is  more  appropriate  for  turbulent  flow. 

3.1.1  Grids 

Tne  grids  used  for  these  calculations  are  rectilinear.  A  portion  of 
the  yrid  for  the  laminar  case  is  shown  in  Figure  3-2.  The  grid  used  for 
turbulent  flow  was  similar,  but  with  finer  spacing  at  the  wall.  Both  grids 
use  geometric  stretching  away  from  the  wall  and  a  constant  spacing  in  the 
streamwise  flow  (x)  direction.  The  acceptable  minimum  y  spacing  (normal  to 
the  wall)  was  determined  by  reducing  the  spacing  at  the  wall  until  the 
solution  was  independent  of  the  wall  spacing.  Details  of  the  grid  are  given 
■in  Tabl  e  3-1 . 


figure  3-2.  Portion  of  rectilinear  grid  used  in  laminar  flow  behind 
a  shock  calculation 
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Figure  3-1.  Schematic  of  boundary  layer  behind  shock:  Region  1  influenced 
by  leading  edge.  Region  2  influenced  by  shock 

the  shock  (Region  2).  As  the  shock  moves  along  the  wall,  both  regions  grow  in 
length  while  the  thickness  of  the  boundary  layer  between  the  two  regions 
i nc  reases . 

One  of  the  distinct  differences  between  the  boundary  layer  immediately 
behind  the  shock  and  the  leading  edge  type  boundary  layer  is  the  direction  of 
the  fluid  velocity  normal  to  the  wall.  In  a  boundary  layer  which  grows  in  the 
direction  of  the  fluid  mean  velocity,  the  fluid  must  flow  away  from  the  wall. 
In  the  boundary  layer  behind  a  shock  the  growth  is  opposite  to  the  direction 
of  the  flow  and  the  fluid  must  move  towards  the  wall  for  the  boundary  layer  to 
grow.  Hence,  the  boundary  layer  edge  normal  velocity  is  negative  near  the 
shock  . 

In  practice,  the  laminar  part  n -  the  boundary  layer  may  be  quite  short. 
However,  the  code  ARC2D  does  not  account  for  transition  from  laminar  to 
turbulent  flow  ( i  .e . ,  computations  are  either  laminar  or  turbulent).  Thus, 


SECTION  3 


TEST  CASES 

Although  the  code  ARC2D  has  Seen  used  for  a  variety  of  flow  problems, 
there  has  not.  been  any  work  with  this  code  on  unsteady  boundary  layer 
development  behind  shocks.  Therefore,  a  study  was  conducted  to  verify  that 
ttie  results  agree  with  experimental  data  and  relevant  analytical  results.  An 
additional  purpose  of  the  study  was  to  evaluate  the  grid  requirements  for 
these  flows.  This  was  done  for  both  laminar  and  turbulent  flows  on  a  flat 
surface.  These  results  are  described  in  Section  3.1.  Calculations  were  also 
made  for  the  case  of  a  shock  interacting  with  a  cylinder,  in  both  inviscid  and 
viscous  flow,  as  summarized  in  Section  3.2.  The  satisfactory  results  from 
these  studies  gives  some  confidence  that  this  method  can  successfully  be  used 
in  appl i cat  ions. 

3.1  BOUNDARY  LAYER  DEVELOPMENT  BEHIND  A  SHOCK  WAVE 

Tne  development  of  the  boundary  layer  behind  a  shock  has  received 
considerable  attention.  Laminar  and  turbulent  flows  have  been  investigated 
noth  analytically4  and  experimentally.^11  The  physical  problem  is  illustrated 
in  figure  3-1.  As  the  shock  wave  moves  along  the  wall  a  boundary  layer 
level  ops  behind  it.  At  the  same  time  the  boundary  layer  is  developing  from 
the  leading  edge  of  the  wall. 

Thus,  the  boundary  layer  has  two  more  or  less  distinct  regions:  that 
influenced  t>y  the  "leading  edge"  (Region  1,  Figure  3-1)  and  that  influenced  by 


viscous  flows.  The  body  surface  pressure  is  calculated  from  the  normal 
momentum  equation,  which  is 


p(l)T’it  +  n x  +  u y )  -  p  U  ( n  x ur  +  OyV£ ) 

=  (nxf, x  +  f,yny)P/r  f  (nx  f  ny)pn  (2-6) 

2  2 
=  Pn  (nx  +  ny) 

where  n  is  the  local  normal  to  the  body  surface.  Free  stream  stagnation 
enthalpy  is  held  constant  along  the  surface  for  inviscid  flows  while  adiabatic 
or  constant  temperature  walls  are  used  for  viscous  flows. 

2.4  TIME  ACCURACY 

Most  of  the  calculations  done  for  this  work  must  be  time-accurate.* 

This  means  time  steps  must  be  limited  to  a  value  that  makes  the  Courant 
number 


r  \  £ 

CN  =  ,  where  c  is  the  velocity  in  the  x-direction, 

less  than  one.  This  requirement,  the  CFL  condition,  means  that  as  the  grid 
becomes  finer,  the  time  step  must  be  reduced.  The  CFL  condition  requires  that 
the  analytic  domain  of  influence  lie  within  the  numerical  domain  of  influence. 
The  maximum  CN  for  all  time-accurate  calculations  for  this  work  is  '-0.9. 


* T n a t  is,  we  are  interested  m  the  transient,  fluid  mechanics,  not  in  obtaining 
a  steady  flow  as  the  time-asymptotic  limit  of  a  transient  flow. 


•  The  region  of  undisturbed  1 1  u  i  <1  that  the  shock  has  not  yet  reached 

•  The  region  behind  the  shock  where  the  fluid  has  been  affected  by 
the  shock  moving  through  it 

'he  computations  described  here  were  done  with  the  following  initial  and 
o  node  ry  cond  i  t  i  ons  . 

Initially  the  flow  is  set  up  with  two  distinct  regimes;  pre-  and 
post  shock ,  as  shown  in  Figure  ?-?.  The  initial  position  and  Mach  number  of 
toe  shock  are  input  to  the  program.  Tha  entire  region  in  front  of  the  shock 
is  assumed  to  no  undisturbed  fluid  with  no  velocity  and  constant  pressure  and 
f  e"iper'  if  ire .  behind  the  shock  the  flow  conditions  are  calculated  using  ideal 
norma  1  shoo:  relationships.  The  inflow  boundary,  7TB,  uses  these  postshock 
'nd’tions  throughout  the  calculation.  If  a  wall  is  adjacent  to  the  inflow 
■ounltry,  the  constant  flow  profile  simulates  physically  a  sharp  leading  edge 
condition.  At  the  "outflow"  boundary,  '(TIT,  the  preshock  conditions  are 
applied,  and  the  calculations  are  stopped  before  the  shock  reaches  this 
noun  i a ry  . 

Along  the  outer  boundary,  FTT,  either  pre-  or  postshock  conditions  are 
up, mied.  Fne  position  of  the  shock  is  calculated  assuming  the  shock  moves  at 
onst ant  velocity  of  Ms  •  aj.  As  one  can  see,  the  actual  position  will  most 
if  fen  :v  between  two  grid  points.  In  the  aljoritnm  used,  the  grid  point 
benini  tne  shore  has  postshock  conditions,  while  the  grid  point  in  front  has 
,>'b‘sri.)cr  tonditions.  Thus,  the  shock  position  on  the  outer  boundary  is 
!>■  f  i  ned  only  to  within  tne  mesh  cell  size.  This  resolution  is  sufficient  for 
n  e <p ! or  if ory  study  such  as  this. 


‘  low. 


different  conditions  at  solid  boundaries  were  applied  depending  on  the 
I  ingeocy  is  satisfied  for  ioviscid  flow,  while  no  slip  is  used  for 


!  he  turbulence  mode  I  in  AkC  XI)  is  .1  conventional  two-layer  algebraic 
n  i  <1  ng  lengt  n  mode  1  .  The  inner  layer  governed  by  the  Prandtl  mixing  length 
w’fn  Van  driest  damping.  The  outer  layer  follows  the  Clauser  approximation. 
Inis  turbulence  model  is  detailed  by  Baldwin  and  Lomax, H  and  is  appropriate 
for  attached  flows.  It.  is  probably  suspect  in  separated  flows. 

►or  AP02D  to  be  employed  in  this  shock -i  nteract  i  on  work,  several 
■ h anges  had  to  oe  made  to  the  code.  The  most  important  of  these  pertained  to 
ranges  m  the  initial  and  boundary  conditions  applied  by  the  code. 

2 .  3  dildNiiARY  AND  INITIAL  CONDITIONS 

One  aspect  of  the  shock-interaction  problem  is  that  the  outer  boundary 
.  mditions  are  time  dependent.  As  shown  in  Figure  2-2,  there  are  two  distinct 
" i ow  regimes  as  the  shock  moves  through  the  computational  domain. 


fi-jure  ?-?.  Schematic  of  initial  conditions  for  blast  wave  flow  problem 


Figure  2-1.  Generalized  curvilinear  coordinate  transformation 

Details  of  this  coordinate  transformation  are  given  in  the  work  of  Pulliam. 3 
The  transformations  are  chosen  so  that  the  grid  spacing  in  computational  space 
is  uniform  with  unit  length.  There  exists  a  one-to-one  correspondence  between 
a  physical  point  in  space  and  a  computational  point  (except  for 
s 1 ngu 1 ari t ies ) .  With  this  construction  a  code  can  handle  a  wide  variety  of 
physical  geometries  with  boundary-conforini ng  grids. 

2.2  NUMERICAL  ALGORITHM 

The  computer  program  used  in  this  work  is  a  version  of  the  code  ARC2D. 
ARC20  was  developed  at  the  NASA  Ames  Research  Center anr|  p, d s  been  applied 
to  a  variety  of  complex  flow  problems.'5**’  ARC2D  uses  an  implicit  finite 
difference  numerical  algorithm  with  fourth  order  dissipation.  The  algorithm 
is  an  implicit  approximate  t  ictori zation  scheme  which  is  second  order  accurate 
in  time.  This  noniterative  scheme  is  based  on  the  work  of  Warming  and  Ream,  ; 
and  results  in  easily  solved  block  tridiagonal  matrices. 

h 
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Fht>  lowest  v  a  I  ijt*  listed,  l .  -I  x  1 1 1  ~  '  ij)pc,irs  to  n>-  due  f  u  a  numerical 
:  illation  in  the  solution  near  the  shock.  Such  osci  I  I  at.  ions  are  present  in 
I  t ne  ;  a  1  mi  I  at i nns  and  are  common  t  i  t  h i  ■ ,  i  inp licit  scheme  with  fourth  order 
a  ■>  i  ■  i  i  ’  i  on  .  Overall  this  study  roved  1  s  the  method  is  well  suited  to 

,  ’  •>..  s'  -  imi.tr  i  ter  of  t  u  rbu  1  e  n  t.  -  v  i  scon  s  effects  behind  the  shock  , 
i*  ••  Cin.c.r  predictions  may  require  a  extremely  fine  grid 

•  ■  -  ■  : :  >q 

•  >w  is  established  when  a  shock  wave  strikes  an  object, 

v  ‘eat  o',.,  h  ts  --Elected  shocks  and  Mach  stems.  To  verify  that  the 
a  a!"  i-  predict  such  flows,  the  interaction  of  a  shock  travelling  past 
ylinde-  was  studied,  both  inviscid  and  turbulent -viscous  flows  were 
1  lilted  to  illustrate  any  differences  in  the  flow  due  to  turbulent-viscous 
f  0  •:  .  *  S  • 

Inv'scid  calculations  of  this  problem  have  been  performed  by  both 
r  1  •  r  and  fernquist2  and  Mark  and  Kutler.15  It  is  important  to  keep  in  mind 
at  tne  inviscid  flow  over  a  ha  1 f-cyl inder  on  a  flat  surface  (the  field 
-jofure  case  for  this  work)  is  identical  to  the  inviscid  flow  over  a  "free 
-inding"  cylinder.  Thus  the  inviscid  shock -cyl  i  nder  flow  represents  both  a 
st  case  and  case  to  compare  with  the  flow  over  a  ha  1 f -cyl i nder  on  a 
rf ice . 

An  analytically  described  mesh  in  a  cylindrical  coordinate  system  was 
ployed  for  the  shock -cyl i nder  interaction  flow  as  shown  in  Figure  3-7.  The 
id  is  determined  by  the  cylinder  body  radius,  the  distance  from  the  body  to 
e  outer  boundary,  and  the  distance  from  the  body  and  the  first  point  away 


urn  the  body,  bays  were  equally  spaced  around  the  body,  and  points  along 
y s  wit"  radially  clustered  near  the  body  by  in  exponential  function.  The 


two  grids  shown  in  Figure  3-/  differ  in  that  the  turbulent  viscous  flow  grid 
has  ,1  larger  number  of  grid  points  along  the  rays  away  from  the  body  with  a 
smaller  distance  from  the  body  to  the  first  grid  point.  The  details  of  the 
grids  are  given  in  Table  3-1. 
i .  2 .  1  Results 

[he  conditions  chosen  for  the  cylinder  flows  were  designed  to  match  the 
experimenta1  conditions  of  Pearson,  et  al.,^  at  the  Ballistic  Research 
Laboratory.  In  their  work  a  0.305m  diameter  aluminum  cylinder  was  rigidly 
fixed  across  a  2.44m  diameter  shock  tube.  A  flattop  shock  wave  with  an 
overpressure  of  42  kPa  (6.09  psi )  passed  over  the  cylinder.  Ambient 
(pres hock )  pressure  for  this  experiment  was  101.3  kPa  (14.69  psi)  with  an 
t mb ient  density  of  1.22  kg/m^.  These  calculations  were  done  with  a  constant 
wall  temperature  for  the  viscous  flow.  The  Reynold  number  based  upon  the 
cylinder  diameter  is  6.95x10^. 

In  many  respects,  the  results  of  the  inviscid  and  turbulent-viscous 
"'ins  ar.;  guite  similar.  Figure  3-8  compares  predicted  pressure-time  histories 
at  several  locations  on  the  cylinder  (u  =  0,  45,  135,  180)  with  those 
experimental ly  measured,  and  with  the  inviscid  calculations  of  Mark  and 
•  itior.  'The  experimental  data  presented  are  an  average  of  the  actual 
measurements  which  contained  substantial  high  frequency  oscillations.)  The 
snarp  rise  in  pressure  at  each  location  signifies  the  arrival  of  the  shock. 

I  be  pressure  at  all  stations,  except  180°,  then  decreases  until  it  "jumps" 
iquin  and  then  reaches  a  plateau  level.  This  second  jump  is  the  result  of  a 
refi<N.?ed  shock  which  is  created  at  the  rear  of  the  cylinder  and  propagates 
i;>st. re,n<i  >Ver  the  cylinder,  as  illustrated  by  the  isopynics  of  Figures  3-9 
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H  jure  3-3.  Pressure-time  histories  for  shock -cyl  i  rider  interaction 
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Figure  3-10.  Computed  turbulent-viscous  isopycnics  of  shock -cyl  i  nder 
interaction  at  various  times 


These  figures  illustrate  some  interesting  features  of  the  computations, 
first  note  that  the  shock  is  inc reasi nyly  "smeared"  far  from  the  cylinder. 

This  is  a  result  of  the  grids.  From  Figure  3-7,  it  is  obvious  that  the 
c i  rvumterent i a  1  grid  spacing  increases  with  distance  from  the  cylinder.  Thus 
the  coarse  grid  results  in  poor  resolution  of  the  shock.  Secondly,  the 
isopycnics  of  the  viscous  flow  are  much  sharper  than  those  of  the  inviscid 
calculation.  This  is  a  result  of  the  finer  mesh  of  the  grid  used  the  viscous 
calculation .  In  addition,  the  increased  damping  of  the  viscous  terms  smooths 
the  isopycnics. 

Noticeable  differences  between  the  experimental  data  and  the 
calculations  in  the  time-pressure  histories  exist  at  the  stagnation  point  (0°) 
and  at  180''.  One  difference  is  that  the  plateau  level  of  the  viscous 
calculation  of  p/pj  -  1.49  is  closer  to  the  steady  state  stagnation  pressure 
ratio  of  1.46  than  is  the  experimental  value  of  1.59.  This  occurs  at  a  time 
when  there  is  a  possibility  of  reflections  from  the  walls  of  the  shock  tube 
affecting  the  measured  pressure.  Therefore,  we  have  some  reservations  about 
the  interpretation  of  the  data. 

Another  possible  source  of  the  discrepancy  at  18(1“  could  be  the 
distortion  of  the  cylinder  during  the  experiment.  The  object  of  the 
experimental  study  was  to  determine  ire  '  .a  •  ■  generi.  aircraft 

structures.  Thus,  the  cylinder  wa1  !  w,-t  •,  *  r  i  >  . :  r  -nun  mi  post  test  inspection 

showed  tnat  the  cylinder  was  deformed .  r.  ,  •  '  .  ;  )  •.  > '  >  that  the 

cylinder  ie  turned  enough  'luring  the  ;  ■  -e.  •*  *  !•  v  *  a  ‘  •  '  ‘ 

‘1  )W*  1  e ’it  . 

f  igur.*  i-Hd  shows  a  growing  t .. -.por.i '•  >m  illation  mi  the  predicted 
i n  v i sc  id  pressure  at  the  IIP  location.  The  sour,  e  of  this  instability  in  the 


inviscid  calculation  is  unknown.  When  viscosity  was  "turned  on"  it 
disappeared  and  therefore  was  not  investigated  further. 

An  important  point  with  regard  to  the  turbulent  viscous  calculation  is 
the  grid  spacing  at  the  cylinder  surface,  fven  though  our  flat  plate  studies 
suggest  a  minimum  spacing  on  the  order  of  10_6,  computation  time  requirements 
precluded  using  a  grid  with  so  fine  a  spacing.  The  minimum  spacing  of  the 
cylinder  grid  (&r/d)  is  HxlO'^ .  With  this  grid,  140,000  steps  are  required 
for  this  flow  calculation  using  11  hr  on  the  CRAY  X-MP  computer.  Thus,  while 
the  primary  features  of  the  flow  appear  correct,  the  magnitude  of  the  friction 
component  of  the  drag  cannot  he  expected  to  be  correct.  Still,  the  calculated 
flow  does  separate  from  the  rear  of  the  cylinder.  Figure  3-11  shows  the 
development  of  the  separation  bubble  behind  the  cylinder.  The  bubble  was 
still  growing  when  the  calculation  was  stopped.  For  comparison  the  stream 
function  contours  of  the  inviscid  calculation  are  shown  .  Figure  3-12.  As 
the  study  of  steady  flows  is  not  the  main  objective  of  this  work,  both  the 
viscous  and  inviscid  calculations  were  stopped  before  steady  state  conditions 
were  obtained. 
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bevel  opnent  of  separation  bubble  behind  cylinder  in  turbulent 
viscous  flow  of  shock-cylinder  interaction  calculation 


SECTION  4 


SHOCK  STRUCTURE  FLOWF 1ELD 


4.1  OH LEM  DESCRIPTION 

Toe  final  study  of  this  work  is  a  calculation  of  the  flow  behind  a 
a  la  st  wave  passing  over  an  object  on  the  earth's  surface.  It  has  a  number  of 
feat  ires  which  are  combinations  of  those  in  the  other  test  cases;  a  shock  wave 
moviny  past  a  cylinder  and  a  boundary  layer  on  a  flat  surface.  We  chose  a  2  m 
radius  half-cylinder  yeometry  to  model  this  problem.  A  consideration  of  some 
other  physical  parameters  reveals  interesting  information,  as  described 
S'  1 ow. 

At  a  distance  of  2000m  from  a  1  megaton  explosion,  the  resulting  blast 
wave  has  approximately  a  1/2  KPa  (25  psi )  overpressure  (Ms  =  1.58).  In  order 
to  numerically  resolve  the  turbulent  boundary  layer  which  forms  behind  the 
ilast  wave,  the  grid  needs  at  least  a  few  points  in  the  viscous  sublayer  (as 
discussed  in  Section  3.1).  At  10m  behind  such  a  blast  wave,  the  turbulent 
wundary  layer  is  less  than  10  cm  high.  Thus,  as  previously  discussed,  grid 
spacing  on  the  order  of  10-&m  is  required.  For  a  calculation  which  is  time 
ac.urjte,  i.e.,  to  satiety  the  CFL  condition,  this  translates  into 
approximately  10 ^  time  steps  for  the  early  time  interaction  with  an  object. 
Such  a  calculation  would  require  over  150  hr  on  the  CRAY  X-MP.  For  this 
present  work,  this  amount  of  computat i ona 1  time  is  unacceptable.  Therefore, 
the  calculation  is  done  with  a  courser  grid  so  that  time  accurate 
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il  illation-,  in.*  made  with  a  reasonable  amount  of  computer  resources. 

'■la soil  on  rue  test  case  results,  the  results  should  he  valid  for  evaluating 
: 1  iha  1  f  I  ow  Mold  e  fleets,  hut.  skin  friction  and  heat  transfer  calculations  are 

not  exported  to  he  !..■  ir  t'  e  . 

Several  quasi -steady  state  calculations  were  attempted  to  supplement 
me  ’■'me  accurate  ca  1 1 1  at  i  on  .  These  runs  included  a  boundary  layer  profile 
t,  rue  inflow  boundary  conditions.  The  thickness  of  these  profiles  was  on  the 
or  lor  of  toe  no-  4h"  >f  t  ue  ,wi  i  ruler,  thus  mode  1  i  ny  conditions  when  the  shock 
,  *  ir  luwusr  re.i.T)  of  toe  protruding  cylinder.  These  calculations  are  not 
i  i  f,  >*d  by  the  small  step  size  requirements  of  time  accurate  calculations,  and 
i •'  f .beret ore  include  t,ne  fine  mesh  required  for  accurate  resolution  of 
/ 1  its  of  fe(;  f  x  . 

'*nt  >r-r .mate  ly ,  these  quasi-steady  calculations  failed  to  converge.  Our 
it’os  wl-  :.ate  that  this  lack  of  convergence  is  due  to  either  the  boundary 
son'  i  *  i  mu  )«-  the  grid  (especially  the  large  aspect  ratio  of  some  cells,  a 
no  ,  ii*  uf  .  i-ip  ifor  limitations).  However,  we  do  not  yet  have  a  definitive 

i  n  s  we  r  '  f  no  goes  t  i  on  . 

\  m  >  lJ  '  ' 

I"-1  jnq  used  in  the  time-accurate  calculations  was  generated  with  a 
nyi'oroo! \  grid  generator  developed  at  NASA  Ames  Research  Center.”  This  grid 
■  >r  produces  a  grid  which  is  orthogonal  to  the  wall,  while  allowing  a 

•  i •  i i  '•>  of  tue  first,  grid  space.  The  stretching  away  from  the  wall  is 
o  ■ >n.oi‘  ill.  The  grid  developed  for  the  time  accurate  calculations  is  shown 

•  •  ;  ,  r .  >  1 .  1  . 


A  •  a  i  eolation  was  performed  t  >r  a  shock  with  a  pressure  ratio  of  ?.  7f> 
■to;:  eg  over  m  inched  s!  mu  i  jr.-  na 1  t  vl  inner)  on  a  wall  (i.e.,  the  ground). 


. . 

. . -f . ! -: . ^0^$^‘:jl':  1 1  • 

■“ivi *Wi'\/\£: >v ■  -..V’1  '■•,[  «\  ■  '•  7T-;i:‘ 

,  \\  -A  _ . JS'-tf . : . . . . rAllr’ 


'  .  . ‘■-tr&A.  '!/'! 

■■'A . . ; . . 


. _ “  f.-V :  ' 


i/.'  I : !  t  Iff; 


:  i  :  :  !  :  i  i  ! 


i  ;,r..  4-1.  i  r  i  >  1  used  for  shock  field  structure  flow 


t  ion  ;.,ot  ioued  until  fne  shock  reached  the  downstream  grid  boundary 
>  ven  cylinder  rdtli  i  towns  f  ream  from  the  cylinder.  A  calculation 
this  was  ,i  ''•formed  b,  Kurylo  ef  al  .  ^  Their  flow  configuration 
f  i  shock  wav-*  w i  *  o  a  1.13  pressure  ratio  passing  an  arched 
’  t  n  a  -it.  i  jv>r.  -to-wi  It-'  ratio  ot  l.  1.  Tne  calculations  of  Kurylo^ 
"led  m  i  -Mar-,,'  ;rid  and  were  stated  to  have  been  carried  out 

,!.  )*•••*."  Howe  v-'  - ,  no  i  nter-ned  i  ate  time  results  were  presented,  and 
COcir  i  j  r  i  t  is  ?.;••)  < :  ■ ; .  i  r  s  e  t  :j  resol  v--  the  viscous  1  yyers  in  the  flow, 
e  tne  difference;  between  their  results  and  ours  be  1  ow. 
re;  i  1 1  s  of  f  he  present  r a  I  c a i at i on  are  presented  in  several  ways. 


1  i  stones  at  different  locations  on  the  cylinder  are  given 


i  he  !  .)<  <i  f.  i  oris  are  idenM  a  1  t  o  those  plotted  tor  the  shock- 
ferae  Mon.  A I  though  'he  shm  k  is  much  stronger  in  tins  case,  the 
f  r  ne  pressu  re- 1 1  me  hi  stories  is  similar  to  the  cylinder  results. 
M’luture  at  various  tidies  is  shown  by  isopycnics  in  Figure  4-3. 
lowtieid  structure  is  similar  to  that  of  the  flow  without  the 
so  nof.>  that  the  shock  does  not  become  more  "smeared"  away  from  the 
result  is  because  the  grid  spacing  does  not  increase  in  the 
;■  i  rvc  1. 1  on  as  i  t  did  for  the  grids  used  in  the  the  shock-cylinder 
A  .11  stimuli  shiny  feature  of  this  flow  is  the  gradual  drop  in 
Mr  is  shown  in  Figure  4-2  M  *  180")  at  t  ~20  ms.  This  drop  in 
r s ,jon  i s  to  the  oev •>  1  opnen t  of  the  separated  flow  as  shown  in 

■ti  the  viscous  Mow  over  the  free-standing  cylinder,  a  region  of 
iw  also  develops  behind  the  half  cylinder.  In  this  case  the 

•  •  waM  accelerates  the  growth  of  the  region.  Here  the 
i*  Mow  *  rom  the  invi  sc  id  solution  is  obvious. 

•  feat  jr.*  which  is  perhaps  more  interesting  accentuates  the 
vi,  ouc  wall  tlow  on  the  flowfield.  Figure  4-5  shows  the 
•  i  s.  guv  uted  Mow  region  in  front  of  the  structure.  This 

•  >•  Mmio.'  "i  the  calculations  of  Kurylo,  perhaps  due  to  the 

• ;  n  *  -M-wdfn  rati)  of  tnei  r  structure.  By  the  time  the  shock  is 
-■Mr s  pas*  the  hal  f-cyl  inder  the  height  of  the  separated  region  is 
e r  iuif.  if  ‘.no  s’r'uct  ire  height  . 

i*i  ir  /,  the  inclusion  o*  v  i  scons- turbi.il  ent  effects  significantly 
*1  iw  ovr  a  hu 1  f  cyUnder  on  a  surface.  The  effects  of  viscosity 
t.  ins,)’  to.  t  he  ne.-g  t^r  d  course  grid,  which  can  not  be  expected  to 
' * > , i 1  •/ . >  the  boundary  layer  near  a  wall  . 
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Figure  4-2. 


Calculated  pressure -time  histories  for  shock-field 
St  rue  tor  interaction 
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a )  14.6  ms 


b)  24.1  ms 


f  i  iju re  4- 


o)  71.1  ms  d )  35. 7  ms 

4.  Stream  function  contours  showing  development  of  separated 
flow  at  rear  of  structure 
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c )  7.1.1  rns 


d)  3b. 7  ms 


Figure  4- 


Stream  function  contours  showing  development  of  separated 
flow  in  front  of  structure 


SECTION  5 


SUMMARY  AND  CONCLUSIONS 

The  flow  resulting  from  the  interaction  of  a  planar  shock  wave  and  an 
object  on  the  ground  was  calculated  with  the  full  Navier-Stokes  equations, 
including  turbulence  effects.  The  resulting  flowfields  show  the  development 
of  regions  of  separation  both  before  and  aft  of  the  structure,  thus 
illustrating  the  importance  of  including  turbulence  and  viscous  effects  in 
predicting  flowfields  resulting  from  blast  waves  passing  over  structures. 

These  results  suffer  from  a  need  to  use  a  coarse  grid  that  can  not  adequately 
resolve  the  details  of  the  boundary  layer  wall  region.  The  use  of  a  finer 
grid  was  not  feasible  due  to  the  inordinately  long  computing  times  which  would 
result,  even  on  the  CRAY  X-MP  computer. 

Test  cases  verified  that  the  implicit  finite  difference  code,  ARC2D, 
could  accurately  predict  boundary  layer  phenomena  and  shock-object  interaction 
given  a  sufficiently  fine  grid.  Still  the  code  does  have  its  limitations  and 
more  effort  is  required  to  solve  the  quasi -steady  flow  over  the  structure  on  a 
wall  (part  of  this  effort  would  probably  be  an  improved  grid-). 

Improvements  in  several  areas  would  greatly  increase  the  accuracy  of 
future  turbul ent-vi scous  flow  calculations.  First,  the  addition  of  some 
scheme  to  the  code  to  reduce  pressure  overshoot  and  the  resulting  oscillations 
would  be  of  benefit.  Second  order  dissipation  near  a  shock  has  been  added  to 
ARC2D  for  time  asymptotic  calculations  and  its  implementation  for  time 


accurate  calculations  should  be  forthcoming.  This  should  be  applied  to  the 
present  cases.  Additionally,  validated  methods  for  calculating  separated 
turbulent  boundary  layers  with  a  course  grid  would  be  helpful  (perhaps 
something  like  a  "wall  function"  for  separated  flow). 
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